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Abstract 

Using the Nambu-Jona-Lasinio model to describe the nucleon as 
a quark-diquark state, we discuss the stabiUty of nuclear matter in 
a hybrid model for the ground state at finite nucleon density. It is 
shown that a simple extension of the model to simulate the effects of 
confinement leads to a scalar polarizability of the nucleon. This, in 
turn, leads to a less attractive effective interaction between the nucle- 
ons, helping to achieve saturation of the nuclear matter ground state. 
It is also pointed out that that the same effect naturally leads to a 
suppression of "Z-graph" contributions with increasing scalar poten- 
tial. 

PACS: 12.39 Fe (chiral Lagrangians), 12.39 Ki (relativistic quark model), 
21.65-|-f (nuclear matter), 24.85-|-p (quarks, gluons and QCD in nuclei 
and nuclear processes). 

Keywords: quark-diquark structure of the nucleon, nuclear matter 
stability, Nambu-Jona-Lasinio model for quarks. 

1 Introduction 

Since the first applications of models based on the linear realization of chi- 
ral symmetry to the description of nuclear many-body systems systems [0], 
the problem of the saturation of the nuclear matter ground state in a ||^ or 
Nambu-Jona-Lasinio 0| (NJL) type models has been extensively discussed 
[i]"0- The common feature of these models is a vacuum effective potential, 
which has a "Mexican hat" shape as a function of the classical scalar field. 
The second derivative (curvature) of this effective potential decreases as one 
moves away from the minimum (vacuum state) towards smaller scalar fields. 
This decrease of the curvature implies an attractive contribution to the effec- 
tive a meson mass at finite density as a result of vacuum fluctuation effects, 
and leads to an increase of the attraction between nucleons associated with 
a meson exchange. These attractive vacuum fluctuation effects, which ap- 
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pear in the form of tadpole type diagrams in a-type models, do not appear 
in non-chiral models like the phenomenologically successful Walecka model 
1^. They lead to an instability of the nuclear matter ground state in the 



mean field (Hartree) approximation. Besides the possibilities of including 
higher order interaction terms in the Lagrangian 0, 0, or fluctuation terms 



arising from the Dirac sea of the nucleons [0, |ri|], no practical methods to- 
wards a solution of this problem have been found so far[]. Since the concept 
of the spontaneous breaking of the chiral symmetry is exhibited most clearly 
and transparently in models based on the linear realization of the symme- 
try, some hints towards a solution of this long standing problem of matter 
stability would be highly welcome. 

In recent years it has been shown that the Faddeev approach to the 
nucleon in the NJL model is a very powerful method to achieve a covariant 
description of the nucleon in terms of quarks and composite diquarks |13|, 
p^ . In particular, this method seems to work very well for the description 
of the form factors and structure functions |]TB[ of a free nucleon. It 
would be interesting to extend this successful description of the free nucleon 
to the investigation of medium modifications of nucleon properties, which 
is currently a very active field both experimentally and theoretically - an 
exciting example are the recent studies of nucleon form factors in (e, e'p) 
reactions []T^. For this purpose, however, one needs the equation of state, and 
one first has to solve the problem of the stability of nuclear matter. A natural 
question in this respect is of course whether or not the quark substructure of 



^Recently it has been argued ||] that the stable abnormal state of quark matter, which 
can be described in the NJL model, should be interpreted as a droplet of massless quarks 
surrounded by the vacuum, i.e., as the nucleon. While this viewpoint is interesting, it 
seems very difficult to proceed to the description of nuclear matter along these lines. 
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the nucleoli plays some role in producing saturation. Important hints that 
this might be so come from the successful quark meson coupling model of 
Guichon and collaborators [jl8[, where the saturation arises from a relativistic 
effect associated with the quark motion inside the nucleon. 

Concerning the energy of a nucleon in a self consistent scalar mean field, 



it has been pointed out [|T9| that there exists an important difference between 
an elementary and a composite nucleon. To second order in the scalar field, 
there is a repulsive contribution proportional to the square of the nucleon 
momentum in both cases. For the elementary nucleon case, this is the famous 
"Z-graph" contribution, but the presence of such a term is actually model 
independent as has been emphasized in a series of papers [0, It is 

very important for the description of nucleon-nucleus scattering and for the 



description of saturation in the Walecka model ||2T| . It is also present in chiral 
models, but at normal densities its contribution to the effective potential 
for nuclear matter is not large enough to stabilize the system against the 
attractive vacuum fluctuation effects discussed above. Rather, the "Z-graph" 
contribution favors the formation of a stable abnormal state since it increases 
strongly with increasing scalar potential (decreasing fermion mass). 

For the case of a composite nucleon, however, there is in general an ad- 
ditional effective, contact-like, a'^N'^ coupling term, corresponding to the 
"scalar polarizability" of the nucleon. If this term is repulsive, it might play 
a role to stabilize the system. To investigate this possibility in a simple NJL 
model calculation is the main motivation for our present work. We demon- 
strate that, as long as there are unphysical thresholds in the NJL model for 
the decay of the quark-diquark bound state, this scalar polarizability is ei- 
ther too small or even of the wrong sign (attractive). However, a sufficiently 
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strong repulsive contribution can arise if confinement effects are incorporated 
such as to avoid the unphysical thresholds. We show that such a term can 
lead to a substantial reduction of the attraction arising from composite a 
meson exchange between two nucleons and therefore to saturation of the 
nuclear matter equation of state. We also demonstrate that the same mech- 
anism naturally leads to a suppression of the "Z-graph" contribution with 
increasing scalar potential (decreasing fermion mass). 

The rest of this paper is organized as follows: In sect. 2 we discuss a hy- 
brid model for the nuclear matter ground state based on the quark-diquark 
picture of a single nucleon. Several expansions in powers of the density are 
given, and the effective two-nucleon interaction (Landau-Migdal interaction) 
is discussed. The numerical results in sect. 3 show that, as long as the model 
contains unphysical thresholds, it leads to essentially the same situation as for 
elementary nucleons; i.e., it cannot describe saturating nuclear matter. By 
using a simple method, based on an infrared cut-off, to avoid the unphysical 
thresholds [^], it will be shown that the scalar polarizability of the nucleon 
discussed above does indeed work towards the stabilization of the system. 
However, for technical reasons, in this work we will limit ourselves to the 
"static approximation" [|^, ^ to the full Faddeev equation when describ- 
ing the nucleon as a quark-diquark state. Since this approximation becomes 
worse for a nucleon in the medium there remains an ambiguity which should 
be resolved in a future numerical analysis. Nevertheless, our results indi- 
cate that a more complete calculation, describing the nucleon in the scalar 
field by the full Faddeev equation, will very probably lead to a saturating 
nuclear matter ground state. Finally, in sect. 4 we present discussions and 
conclusions. 
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2 Nuclear matter in the NJL model 



The NJL model is characterized by a chirally symmetric 4-fermi interaction 
between quarks. Any such 4-fermi interaction can be Fierz symmetrized [|12[ 



and rewritten identically as a chirally symmetric linear combination of the 
to™ E. O. (?^.,)^ whe. , is the flavo. SU,2) ^ua.. «e.d. P. a. ,„at*es 

in Dirac, flavor and color space, and the coupling constants, Gi, are functions 
of the coupling constants appearing in the original interaction Lagrangian. 
Writing out explicitly only the scalar, pseudoscalar and vector terms, which 
are relevant for our present discussion, we have 

£ = ^(^^ - m) + ((V^^)' - {Hl5r)^y) - (i^Ytpf + ... (2.1) 

where m is the current quark mass. In the nuclear medium characterized 
by the density p, the quark bilinears tpip and ipj^tp have expectation values 
which we will separate as usual according to i/ji/j = (plipiplp) + (: ipip '■) 
ip'^t^ip = {plip'y'^iplp) + (: ■ip'^^^ip :). The Lagrangian can then be expressed as 

C = ^{t^-M- 2G^Yu^,) ^ - ^ ' + G^u^u^ + (2.2) 

where we defined M = m — 2G.„{p\il)il)\p) and cu^ = (p|'?/'7'^'?/'|p), and Ci is 
the normal ordered interaction Lagrangian^. 

•^Our conventions are as follows: In the effective potential, M is considered as a variable, 
and its 'physical' value is denoted as M* (solution of the gap equation in the medium) 
or Mq (solution of the gap equation in the vacuum). The quantity {p\i}}ip\p) is considered 
as a function of M, and its value at density p and M = M* is denoted as (tpip)* , while 
its value for zero density and M — is the physical vacuum value {ijj'4')o- Concerning 
the diquark and nucleon masses Md = Md{M) and Mm = Mn{M), their values at 
M = M* (M = Mo) are denoted as and {Mdo and Mno)- Similar notation 
is used later for the energy density, £, the nucleon energy, e{k), the forward scattering 
amplitude, f{k', k), the TriV sigma term, EttAt, etc. (E.g; £* = £{M = M*), etc.) Further, 
all four dimensional momentum integrals, like that in Eq. (|2.5| ), stand symbolically for their 
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We briefly recapitulate the procedure to construct the nucleon as a quark- 
diquark state for zero baryon density. Using a further Fierz transformation, 
one can decompose the interaction Lagrangian in Eq. (|2.2|) into a sum of qq 
channel interaction terms |]T^. For our purposes we need the term which 
describes the qq interaction in the scalar diquark (J'^ = 0'^,T = 0, color 3) 
channel: 

= Gs (75C) T2P^i^^) (V^^ (C-S5) r2P^i^) , (2.3) 



where (3 = y 3/2 A {A = 2, 5, 7) are the color 3 matrices, and C = 27270- 
The coupling constant Gs is a function of the coupling constants appearing 
in the original interaction Lagrangian. The reduced t-matrix in the scalar 
diquark channel is then obtained from the Bethe-Salpeter (BS) equation as 

1 + 2GsUs{q) 



with the scalar qq bubble graph 



r (Tk 



^,S{k)^,S{-{q-k)) 



(2.5) 



Here 



is the Feynman propagator for the constituent quark, and tr/j denotes the 
trace over the Dirac indices. If one restricts the interacting qq channels to 
the scalar one, the relativistic Faddeev equation [|1^ can be reduced to an 
effective BS equation for a composite scalar diquark and a quark interacting 



via quark exchange ||25|. This has been solved numerically by using the 

regularized expressions in one of the schemes considered in this paper (EucUdea n cu t-off, 
3- momentum cut-off, or proper time cut-off scheme). We also note that £7 in (2^) also 
contains the counter terms linear in the normal ordered products. 
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T{P) = TJ^-^TT^ (2-7) 



Euclidean sharp cut-off scheme. For finite density, however, the gap equation 
(to be given below) depends on the nucleon mass, and in order to find its 
solution one has to solve the Faddeev equation many times. In this paper, as 
a first step, we will restrict ourselves to the static approximation [^, ^ to 
the Faddeev equation, where the Feynman propagator in the quark exchange 
kernel is simply replaced by — 1/M. Then the Faddeev equation reduces to 
a geometric series of quark-diquark bubble graphs (IIn{p)), and the solution 
for the t-matrix in the color singlet channel is 

3 1 

Mi + -|njv(p) 

with 

Il^{p) = -J-0^,S{k)r,{p-k). (2.8) 

The nucleon mass M^v is then obtained as a solution of the equation 
l + -^njv(^=Mjv) =0. 

2.1 A hybrid model for the nuclear matter ground 
state 

Using hadronization and path integral techniques, any quark Lagrangian of 
the NJL type can be expressed in terms of physical hadron degrees of freedom 
(nucleons and mesons), see Ref.[|^, |T^. The result is a complicated non-local 
Lagrangian of the linear a model type, which also contains interaction terms 
of higher order than in the standard a model. Since in the course of the 
derivation one has to integrate over the quark fields as well as the auxiliary 
diquark and color octet baryon fields, "trace log terms" also appear in the 
Lagrangian and these should be taken into account in the computation of 
the effective potential. The resulting a model type Lagrangian can then, in 



principle, be applied to nuclear matter by introducing a chemical potential 
for the nucleons and using some standard approximation scheme like the 
mean field approximation. 

In this paper, however, we will not fully exploit this ambitious method, 
but we will resort to the most simple, nontrivial approximation. First, in the 
"trace log terms" we take into account only the quark loop term, that is we 
neglect the diquark loop and baryon octet loop terms which can be thought 
of as ring type quark-quark and quark-diquark correlations in the Dirac sea 
of quarks. Second, the resulting effective o model-type Lagrangian is treated 
on the mean field level. In this way one obtains the following expression for 
the effective potential (energy density )[]: 

E = Ey- G^uj^uj^ + 77V / ^n(fc)e(fc), (2.9) 

where the contribution due to the quark loop is 

with 7q = 12. In Eq.( |2.9| ), n{k) is the Fermi distribution function for the 
nucleons, = 4, and e{k) is the energy of the nucleon (pole of the quark- 
diquark t-matrix) moving in the scalar and vector fields which couple to the 
quarks. 

This simple model can also be motivated by defining a "hybrid approxi- 
mation" for the nuclear matter ground state as follows. The nuclear matter 
expectation value of any local quark operator, O, consists of its expectation 
value in the "valence quark vacuum", |0) = |p = >, and an average over 

■^In this paper we will frequently refer to the function £ as the "effective potential", 
even though we use the density as a variable rather than the chemical potential. 
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the nucleoli Fermi sea consisting of correlated valence (f ) quarks: 



/fi fC 
j^^n{k){OUk), (2.11) 

where 

{0),{k) = J d^r[{N,k\0{r)\N,k) - (O|0(r)|O)] (2.12) 

with |A^, k) denoting the correlated valence quark state (quark-diquark state 
in our case) of the nucleon. (We divided throughout by an overall volume 
factor V in Eq.( |2.11| ).) Taking O = Ti, where Ti. is the quark Hamiltonian 
density corresponding to the Lagrangian density (|2.2D, we have {H)v{k) = 
e{k), and (0|7i|0) in the mean field approximation is given by Sy — G^uj^oj^ 
of Eq.( |2.9| ). In this way one arrives at Eq. ( |2.9| ). We will see later that, 
for O = ipip or Ti, Eq.( p.ll| ) is consistent with the well known low density 
expansions of the quark condensate and the energy density. 

2.1.1 Mean vector field 

Let us first discuss the form of e{k) and determine the mean vector field: As is 



clear from Eq.( p.2|) , the constituent quark propagator in the scalar and vector 
fields is obtained by replacing k^ kq^ k'^ — 2GujUJ^ in the propagator ( |2.6|) . 
Therefore, the qq bubble graph in the mean field approximation is obtained 
from ( p.5|) by replacing S{k) S^kq) and S{—{q — k)) S{—{q — k)Q) = 
S{kQ — qd), where q*^ = — AG^u'^. Thus the bubble graph and the 
reduced t-matrix in the mean fields are obtained by replacing q'^ in 

( |2.5|) and ( p.4|) , respectively. Consequently, the quark-diquark bubble graph 
in the mean fields is obtained from ( |2.^ ) by replacing S{k) S{kQ) and 
t{p -k) ^ t{{p - k)n) = t{pn - ko), where 

p^=p^-6G^^^ (2.13) 
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As a result, the quark- diquark bubble graph and the quark-diquark t-matrix 
in the mean fields is obtained by replacing p% in ( p.8|) and (|2.7|) , 

respectively. The spectrum is then obtained by noting that if the equation 
1 + j^Il]\r{p) = is satisfied at ^ = Mjy, the solution of 1 + j^II^{pn) = is 
•^N = Misf, and solving for p° = e(p) this gives the positive energy pole as 



e{p) = GG^oj" + Jp% + Ml = 66-^^° + En{pn). (2.14) 



Using this in Eg. ( p.9|) we obtain the physical vector field (u;*'^) from the 
condition {d£ j du^) = as 

UJ*^ = 3j^, (2.15) 

where the baryon current is given by 

■0 • [ 



Using this to eliminate the vector field, we obtain for the energy density 

£ = £v + £u. + £f, (2.17) 

where the vacuum part has been given in Eq. ( |2.10| ), and 

£u = 9G'^(p'+j|), (2.18) 
j^n{k)E^{kM). (2.19) 

Unless stated explicitly, we will consider nuclear matter at rest, where uj* = 
j Q = 0, k]\r = k, and n{k) becomes the usual step function Q{pf — k). 
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2.1.2 Mean scalar field 



The mean scalar field in the medium (or the constituent quark mass) is 
determined by the condition 

M — fn 

dS/dM = + (pI^^Ip) = 0, (2.20) 

where the second equality is the gap equation which determines the physical 
value of M (M = M*). The quark condensate in ( p.2CI| ) is given by (see 
Eq. (CT) ) 

^n(fc)(^^),(fc) (2.21) 

with 

(0|«.|0> = -2n,M/|^j,-^ (2.22) 

) ^^^^^ \ ) 

We finally note the following relation, which follows from ( |2.20|) : 

2^1^ = (2 24) 

'^dm 1 + 2G.(^(0|^V^|0))^ 



M=M* 



(2.25) 



V^^^-^ / M=M* ^^^CT 

Here g is the meson-quark coupling constant at zero density and zero meson 
momentum^, and Eq. ( p.25|) defines the effective a mass in the medium, M* = 
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The meson-quark coupling constant g is introduced such that, in terms of the nor- 
mahzed meson fields a = —^^ipi^ ^^'^ = ~^y^V'^75''"'/'j the Yukawa coupling terms in 
the semi-bosonized Lagrangian read —ipg (cr + ztt • T75) ip. In terms of the pionic bubble 
graph in the vacuum (n^y(/c^)) it is given by c/^ — -i/^§gr- at — and M = Mq. We 
also note that for to = we have the relations Mog = f.^ and M^q = 2Mo from the gap 
equation in the vacuum and the qq bubble graph for pion decay. 
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Mfj{M = M*), in terms of the curvature of the effective potential. The 
zero density analogs to (|2.24 ), ( 2.25 ) are obtained by replacing M* — > Mq, 



£ ^ £v, Ml ^ Ml,. 

2.2 Expansions in powers of the density 

In this subsection we derive various model independent expansions in powers 
of the density, as well as some relations which are peculiar to the hybrid 
model introduced in the previous subsection. These discussions will help to 
understand the physics behind the hybrid model expressions ( |2.9| ) or ( |2.11| ), 



in particular the role of the vacuum terms. They will also be useful when 
we discuss the numerical results in sect. 3, in particular the mechanism for 
saturation. 

The validity of the Feynman-Hellman theorem 

{^^l^Y - (^^)o = (2.26) 

in the present hybrid model is easily seen as follows. Relation ( p.2(]| ) at the 
"physical value", M = M*, reads M* = m — 26"^ (■?/;■?/;)* (see the first footnote 
in this section for our notation). From this we can express the difference 
between the condensate in the medium and in the vacuum as (tpip)* — (iptp) o = 
— (M* — Mq)/2Gt^. On the other hand, since the energy density is stationary 
at M = M* we also have the relation dS*/dm = dS*/dm = — (M* — 
Mq)/2Gt,, where in the last equality we used the form (|2.9|). This leads to 
the Feynman-Hellman theorem ( p.26|) . 

The expansion of the energy density S* = £{M = M*) around some fixed 
value of the density can be derived in a general (model independent) way by 
applying Landau's expression for 68* to the case of spin-isospin independent 
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variations of the nucleoli distribution functions, 6nk, caused by a deviation 
of the Fermi momentum from some value ppo (density Pq): 

(2.27) 

The Landau-Migdal interaction refers to the case \k'\ = \k\ = pf and is 
denoted as /*(cosG), where is the angle between k' and k. We insert 

^"^-t^^^^i-^r-'^Hp.-k) (2.28) 

into Eq.( p.27| ). First we derive the low density expansion and set ppo = 0. 
Note that each phase space factor in (|2.27|) contains a factor k"^, so that the 
first term of (|2.27|) involves p^ and the second term p^. The first term is 
simply the expansion of the energy density of a free Fermi gas (Sq) in powers 
of pp. Performing partial integrations for the second term we obtain the low 
density expansion 

where the Landau-Migdal parameter /£^q is defined as /£^q = ^ j\ d cos9 
/*(cos9), and /o,l=o is the zero density value |. 

^The relation to the spin averaged s-wave scattering length in the vacuum (aL=o) is 
/o,L=o = jv/^p o,L=o- In order to obtain the low density expansion, we assumed the Landau- 
Migdal interaction to be finite. The contribution of the pionic Fock term is proportional 
to the average over the Landau angle of q^/(<j'^ + -^^^o)^ with — 2p|,(l — cos0). If 
one sets M^^q = from the outset, this goes as 1/pp, and its contribution to the energy 
density and the quark condensate would be oc p'^^^ instead of ||2^. However, in the 
real world the low densit y exp ansion always implies pp << M^^o- We also note that the 
origin of the term in ( |2.29| ) is, besides the relativistic correction to the free Fer mi gas 
energy density, the three-body term which is the next term in the expansion ( ^.27| ). The 
three-body reducible piece of this term behaves as l/pj? as pp 0, and therefore gives 
rise to a contribution cx aj^^Q p'p p8| . 
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The corresponding low density expansion of the quark condensate is then 
obtained from the Feynman-Hellman theorem ( 2.26|) as 



dM^o (, 3 pi_\ , p'dfo,L 



<«• - <«» - ^ - lost r TisT + °(^»- (^-^"i 



We introduce the following definitions of the one and two body vrA^ "a-terms" , 
which we give here for the case of finite density for later use0: 

E:-..™^; S:..-"lf. P.31) 

where e*p is the Fermi energy and Vp = {de*{k)/dk)f,^p^ is the Fermi velocity. 
For zero density, the first relation in ( p.31| ) reduces to the usual definition of 
the one body nN cr-term (Sttat^o = ^ ^ dr» ° ) • '^^^ Gell-Mann-Oakes-Renner 
(GOR) relation in the NJL model reads m m (V'V')o = -{F^M^q)'^ C with 
C = (l — where we define the pion decay constant, F^, and pion mass 
in the vacuum, M^ro, at zero momentum instead of at the pion pole. Then 
Eq. ( |2.30|) can be written in the following form: 

~ ^ {F^M^ofC y lOM^o/ 2 {F^M^ofC^^^^^^- ^'''''^ 

All relations given above are model independent. It is a simple task to 
check their validity in our hybrid model, and to give the corresponding low 
density expansions of the constituent quark mass and the effective nucleon 
mass: Starting from the expression ( |2.17| ) for the energy density and using 



^Our notation in (2.31) indicates that the quantity S^*jv is the derivative of the mo- 
mentum dependent one-body cr-term at finite density — m ^^^j;^) at k = pp. We 
remark that, to our knowledge, the nNN sigma term of Eq. ( p.31| ) is introduced here for 
the first time. It is the expectation value of ip^p between a correlated two-nucleon state, 
and it would be interesting to investigate whether there is a relation to observables of 
pion-deuteron scattering. 
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the fact that the vacuum part Ey is stationary at M = Mq, we obtain the 
expansion 



(M - Mq)^ 
AG ^ 

^ dm 



p{M- Mo)^^ + 9G^p' + ..., (2.33) 



where we used ( ^.241 ), ( f^.25D in the vacuum. Minimizing ( f^.33[ ) with respect 



to M we obtain M* — Mq = — 2G'^p ^^" + . . ., which is the lowest order term 
in ( ^.3U| ), and inserting this into the above expression gives 



S*-Sn 



1 2 (OMn^ 

2^ ['m 



+ 9G^p^ + ... (2.34) 



M=Mo 



2P 



9lo 
Ml, 



M2 



+ . . . 



(2.35) 



which is Eg. (|2.29|) with the scattering length given by a and u exchange 
between the nucleons. The factor | arises here because the vacuum part 



(first term in ( |2.33| )) cancels half of the attractive Fermi part (second term 
in ( |2.33| )). The a meson mass has already been defined in terms of the 
curvature of the effective potential in (|2.25| ), and Qao is the vacuum (zero 
density) value of the aNN coupling constant at zero momentum, g*, which 
is defined by 

(2.36) 



N 



DM 



M=M* 



9a 

9 ' 



(The aqq coupling constant at zero density, g, has been introduced in sect. 
2.1.2.) For the lo meson exchange part in ( |2.35[ ), we identified 

1 9l 



Gu 



(2.37) 



18 

Since we have already confirmed the validity of the Feynman-Hellman 
theorem ( p.26D in our hybrid model, the relations (|2.30|) and (|2.32|) are also 



16 



validQ The corresponding expansions of the constituent quark mass and the 
effective nucleon mass in our model are then easily obtained as 



M* 
Mo 



-"n-N,0 



NNfi 



2 {F^M^o? 



M 



N 



M 



NO 



A, 



(F^M^or 

^ttNN,0 



(F^M^o) 




+ Oip 



(2.38) 



(2.39) 



where in (|2.39|) we introduced the zero density values of the following quan- 
tities: 



A* 



dM 



Mo 
Ma 



B* 



(9M2 



Ml 

Mr, 



(2.40) 



The "naive" (Mjv = 3M) values are Ao = A* = 1 and Bq = B* = 0, and for 
this case the term proportional to p in ( ^.391 ) agrees with the result given by 
Eq.(5.16) of Ref.[|0|. 

In a similar way one can expand the energy density, quark condensate, 
etc, around any density (po > 0). The difference is that such an expansion 
in powers of (p — po) is analytic, while the low density expansion discussed 
above is analytic in the Fermi momentum rather than the density. Therefore, 
in order the get the expansion of £* up to second order in (p — po), it is 



^It is, of course, also possible to derive ( 2.3C ) directly from the hybri d model relation 



( 2.21 ). For example, the term oc p is easily obtained by expanding ( 2.21 ) as follows: 



By using M* 
term in ( 2.3C ) 



Mo) 



a(o|vV'|o) 



dM 



P 



N 



M=Ma 



\ dM 



M=Mo 



(V'V)o + (M* 

— Mo = — 2G7r(('0V')* ~ (^V')o) and Eq. (|2.24| ) we easily obtain the linear 
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sufficient to consider the n = 1 and n = 2 terms in ( p.28| ) to obtain the well 
know expression 

r2 \ 



1^ ,2/^' 



S*{p)-S*{po) = (p-po)eF(Po) + 2(p-Po) I 2p|^ 



vUPo) + fl=oiPo)j+0{{p-pof 

(2.41) 



Using the Feynman-Hellman theorem ( ^.26| ), this can be converted into the 
corresponding expansions of the quark condensate, and in our hybrid model 
this in turn gives the expansions of the constituent quark mass and the 
nucleon mass around any density Pq: 



TT 



2p2 



(2.42) 



FO 



M*{p)~M*{po) -1 f, ( ^^^( ^^ 



7r^^7rV(Po) + ^InnIpo) 



(2.43) 



mUp)-mUpo) 



M, 



m 



(p - po)A*(po)s;^(po) + T^(p - Pof 



A*{po 



TX 



Mo 



j:':N{Po)+KNNiPo)]-B*{po 



(2.44) 



We finally note for later use that in our model the one-body cr-term at finite 
density can be expressed as (see eqs. (|2.14|) , (|2.24| ), (|2.25| ) and (|2.36|) ): 



del 



-'ttN 



m 



m- 



N 



dM* 



m gg* 



dm "~E*p\dM J dm 2G^ E*p Mf 



(2.45) 



and therefore the derivative of the constituent quark mass w.r.t. the density 
is determined by the a meson exchange interaction between a quark and 



a nucleoli, while the derivative of the effective nucleoli mass is determined 
by the nucleon-nucleon interaction due to a meson exchange. (Note the 
additional factor A* (x g^/ g in (|2.44|) .) This point will be discussed further 



in the next subsection. 

2.3 The Landau-Migdal interaction and a mechanism 
for saturation in chiral Hartree-type theories 

In this section we derive the form of the effective two-nucleon interaction 
(forward scattering amplitude) from the definition ( |2.27| ), and discuss some 
important features which arise from the quark substructure of the nucleon. 
For this purpose, we follow closely the method explained in Ref. [0. Since 
the energy density is stationary at X = X*, where X denotes both M and 
ui^, the quasiparticle energy is given byQ 

. 5£* dS* 



(2.46) 



* 6ni drii ' 
and the effective two-nucleon interaction by 
_ 6^8* _ 6e* _ de* ^ / geA 



Here we use the relations 

S fdS\ ^ _ / d^£ \ ^ f d^£ ^ dM* 



Srij \dX J \dn.jdX ) ^^-^^ \dXdM J ^^^^ dn 



■J 

[dXdujJ^^^, drij ■ 

In the limit of zero baryon current (nuclear matter at rest) we have [0 
gg§^ ^ 0, and we can use (|g) to eliminate ^ and ^ from Eq.(|^). 



^Here the label i denotes a quasiparticle state. We use the symbol 6 for total variations, 
and d for partial variations. 
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This gives the interaction in the form 



ft 



de* 



drij 



de. 



duJu. 



x=x* 
(2.49) 

where the first term is zero in our hybrid model (see ( |2.14| ), and the other two 
terms have aheady the form of the direct a and u) exchange potentials with 
the meson masses given by the second derivatives of the energy density and 
the couplings to the nucleon by the derivatives of the quasiparticle energy 
with respect to the mean fields (see eqs. (|2.25| ), (|2.36|) for the a meson). 
Carrying out the required derivatives by using ( ^.9] ) and ( |2.14| ) and setting 
= at the end, we obtain in our model for two nucleons with momenta 
p, p' : 

np',p)- 



9, 



*2 



+ 



9. 



9: 



p p 



(2.50) 

where Ep = E^{pp). The aNN coupling constant is given by ( |2.36| ), and 
the squared a mass has been split into a vacuum part and a density part 
{Ml = Mly + IIctd), which are given by 



Mly 



1 



2^7, 



n 



aD 



9llN 



d'^k k^ + M2 



(27r)4 (p - M2 + 



le 



(2.51) 



d^k Mn 
^n{k) 



(27r 



E^ik)- 



(2.52) 



The (TcrNN coupling constant in Eq. (|2.52|) arises due to the curvature of the 
function Mn{M) and is defined as 



9ao_ 

9' 



N 



(2.53) 
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Concerning the uj meson exchange part in ( |2.5CI| ), we used the definition 
( 2.37 ). The last term in ( p. 50 ) contributes to the Landau-Migdal parameter 
fi^ii and the first two terms give 



'M 



N 



L=0 



El 



^ *2 2 



(2.54) 



This resuh, which we derived here for clarity directly from the definition 
of the Landau-Migdal interaction, can of course also be obtained from the 
developments of the previous subsection. From eqs. (|2.43|) , ( p.44|) and ( p.45| ) 
we see that the derivative of the effective quark and nucleon masses w.r.t. 
density can be expressed as 



DM* 
dp 

dp 



El 



'dM, 



N 



M 



N 



E*p 



dM 



[dAP J 



M=M* 



M=M* 



(2.55) 
(2.56) 



and therefore the derivative of the nucleon energy ( p.l4|) at the Fermi surface 

/de*(p)\ 



IS 



dp 



L=05 



(2.57) 



P=PF 



where /£^q is given by (|2.54|) . 

From the developments in this and the previous subsection we can make 
the following important observation. If the function M]\f{M) for the nucleon 
mass in the scalar field has an appreciable positive curvature d^M^/dM'^ for 
large scalar potentials^ (small M), (i.e., if Mn{M) becomes "fiat" for small 
M), there arise two important repulsive effects. First, the aNN coupling 
^Here and in the following, the "scalar potential" is defined as Mq — M. 
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constant ( p.36| ) is reduced for large scalar potentials (large densities). Sec- 
ond, the sigma mass in the medium receives a repulsive contribution from 
the second term in (|2.52|) , which can be expected to over-compensate the 
reduction of the first term at normal densities, since the first term arises 
entirely from the Fermi motion and is of the order 1/M%. Therefore, an 
appreciable curvature of the function M]^{M) could lead to a suppression of 
the attraction due to the a meson exchange contribution to /2=0' (|2.54|) . 
This could lead to a scalar potential for the nucleons which grows less rapidly 
with increasing density than the vector potential (which grows oc p). If the 
r.h.s. of Eq. (|2.56|) is a decreasing function of the density, the scalar poten- 
tial acting on the nucleons (oc M^q — M^) will grow less rapidly than oc p. 
This "decoupling mechanism" of the scalar and vector potentials for higher 
densities is a necessary condition for saturation in any Hartree-type theory 



MM 



To appreciate this point, let us discuss the behavior of the a meson mass 
as a function of the scalar mean field in more detail. As we discussed in sect. 
1, the Mexican hat shape of the vacuum effective potential {Sy) implies that 
the vacuum fluctuations contribute attractively to the sigma mass. This 
can be seen from (|2.51| ), which decreases with decreasing M and becomes 
negative below some value of M (tachyon pole). This leads to an attractive 
contribution to the interaction ( |2.54D , increasing with density. The density 
dependent nucleon loop contribution to the a mass in Eq. (|2.52|) consists of 
two terms. The first one is the Fermi average over the "Z-graph" contribution 
discussed in sect.l. It is repulsive and for the case of elementary nucleons, 
where is a constant independent of the scalar field, increases rapidly with 
increasing scalar potential (decreasing M^r). Nevertheless, from a and NJL 
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model studies for elementary fermions it is known that, irrespective 

of the model parameters, the size of this term is not sufficient to stabilize 
the state at normal densities - although for higher densities it does tend to 
stabilize the abnormal state. 



The second term in ( |2.52| ) is the Fermi average over a auNN contact- 



type interaction, which arises from the "scalar polarizability" (oc ^g^i ) of 
the nucleon - i.e., the response of the internal structure of the nucleon to the 
applied scalar field. The presence of such a term, which vanishes in the case of 



elementary fermions, has been pointed out in a series of papers [|19i pO| . If the 
curvature of the function Mjv(M) is positive everywhere and increasing with 
decreasing M, we can expect that this term gives a stabilizing contribution 
which might lead to a stable normal state. In this connection it should be 
emphasized that confinement effects are expected to lead to a nucleon mass 
which does not vanish as the quark mass M vanishes. That is, confinement 
effects should make the curve Mn{M) "flatter" than the naive case = 
3M, in particular for small M. This would imply a positive curvature which 
increases, and a slope which decreases, with increasing scalar potential - i.e., 
a repulsive contribution due to the contact term and a suppression of the "Z- 
graph" . This important point will be investigated further in a simple model 
calculation in the next section. 

2.4 Chiral symmetry and connection to the linear sigma 
model 

The contributions to the a mass in the medium can be expressed in terms 
of the derivative of the quark condensate with respect to M, see Eq. (p.25|) . 
For the pion, the analogous relation is obtained from the low energy theorem 
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{p\ipip\p) = MIlT^{k = 0), which holds in the NJL model between the quark 
condensate and the pionic polarization operator at zero momentum (see e.g; 



Appendix B of Ref. ||29||). Together with Eq.( p.20D , this shows that the pion 
mass in the medium, = g'^ (^2k~ ~^ ^iri^ = 0)) , is related to the derivative 
of the energy density as 

d£ _MMl m 

DM ~ 2G;- ^^-^^^ 

At the physical point M = M* and for exact chiral symmetry (m = 0), one 
therefore has the possibilities M* = (Goldstone mode) or M* = (Wigner 
mode). Explicitly we have from eqs. (|2.21| )-( p.23| ): 

I" d?k 1 
+ (sJ + 9..)7./p^«(*)^. (2.60) 

where ( p. 591) is the vacuum term (M^y), which leads to a tachyonic pion for 
M < Mq, and for the Fermi sea contributions we introduced the nNN and 



TTTcNN coupling constants 



9 M 
g^^ Mm (OMm M^ 



(2.61) 
(2.62) 



^2 M \dM M )' 

The contribution proportional to g1 in ( ^.6(J| ) is the Fermi average over the 
"Z-graph" contribution with external zero momentum pion linesQ, and the 

^"This follows from the fact that the chiral Ward identity implies r7r(0) = gjbT for 
the ttNN vertex at zero pion momentum. We also note that the coupling constant g^^ is not 
the physical ttNN coupling constant, but it corresponds to the coupling constant appearing 
in the tt NN and a NN interaction Lagrangians of the linear a model for elementary 
nucleons |Q. In the linear a model for nucleons, the coupling constant appearing in the 
Lagrangian (denoted as here) is defined [ |3l[ | by the relation Mjvo = where v is the 
vacuum expectation value of the renormalized a field, which in our notation is given by 
V = M/g. 
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term proportional to g-jm, which vanishes in the "naive" case 3M, is 

the Fermi average over the vrvrNN contact term. It should be noted that the 
pionic "Z-graph" is not suppressed (but rather enhanced) by the mechanism 
discussed in the previous subsection, but nevertheless the sum of the Z-graph 
and the contact term, which is proportional to might be suppressed 

if the curve Mjv(M) becomes "flat" for small M. 

Using the expressions (|2.25| ) and ( p.58|) for the derivatives of the energy 



density in terms of the sigma and pion masses, it is easy to show that that 
vacuum part, Sy, in Eq.( p.lO ) is essentially equivalent to the "Mexican hat" 



of the linear sigma model. Using ( |2.10| ) we can write Sv{M) = Si{M'^) — {M— 
Mo) where Si depends only on M^. Expanding Si around = Mq we 
obtain 



SviM) = -{M 



M2=M2 



Since the derivatives of Si at M = Mq are expressed as 



(2.63) 

'dsA _ mmIq 



(^\ _ Ml, 



and TTTTTT = — TTi the above expression becomes 

\dM^)M=M, 9 



m /_o , .ox2 Ml - M3n /, ,o , Ml^ 

-,2 



Sy{M) = - (M - Mo) ^ + (M^ - MlY + (M^ - M^, ^ 



+ O i^[M' -M^) j . (2.64) 

Except for the terms of third and higher order in (M^ — Mq), this is the 
Mexican hat of the linear sigma model. 
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3 Numerical results and discussions 



Let us now discuss the numerical results based on the expression ( 2.17 ) for 
the energy density. First we refer to a calculation where confinement effects 
(removal of unphysical thresholds) are not incorporated. The parameters of 
the model are the four-fermi coupling constants G.,,, Gg and G^, the cur- 
rent quark mass m, and the ultra violet (UV) cut-off A. Gtt, m and A are 
determined as usual [0 by the conditions m^=140 MeV, fn=d3 MeV, and 
Mo=400 MeV. (Here m^r and f-,, are defined at the pion pole, in contrast to 
M^o and used earlier. We also note that to fix the solution of the gap 
equation at zero density as Mo=400 MeV is rather arbitrary, but qualita- 
tively similar results are obtained for other reasonable choices.) Lastly, Gg 
(or the ratio Vs = Gs/G^) is determined by requiring MAro=939 MeV for the 
pole position of the quark-diquark t-matrix in the static approximation (|2.7] ) 
at zero density. In particular, for all cases considered here this requirement 
implies a pole of the diquark t-matrix ( |2.4| ) at the scalar di quark mass M^, 
and the numerical calculation shows that this pole term, together with the el- 
ementary 4-fermi interaction {4iGs), overwhelms the contribution due to the 
qq continuum. We will therefore, for simplicity, approximate the qq t-matrix 
a-s Ts(q') — 4:iGs + gs/{<f ~ ^d) following calculations, where Qs is the 

residue of ( |2.4| ) at the pole [Q. The parameters determined in this way for 
the Euclidean sharp cut-off, the 3-momentum sharp cut-off, and the proper- 
time cut-off schemes are listed in Table 1. G^ is the only parameter which 
is left for the finite density calculation, and for all cases considered here we 
fix it by the requirement that the curve for the binding energy per nucleon 
{Eb/A) passes through the point {p,Eb/A) = (0.16/m-^ -15MeV^). We 
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first discuss the case of tlie Euclidean sliarp cut-off, since for this case the 
exact Faddeev results for the nucleon mass are available [0 for comparison. 

Fig.l shows the diquark mass (dot-dashed line), the quark-diquark thresh- 
old M/) + M (upper double-dot dashed line), and the nucleon mass in the 
static approximation (lower dashed line) as functions of the scalar potential 
$ = Mo - M. (The dotted line shows M = Mq - $ for orientation.) For 
comparison we also show the result of the exact Faddeev calculation for the 
nucleon mass [0, ^ by the solid line. (In the exact Faddeev calculation, 
a slightly larger ratio, = 0.655, has to be used in order to get the same 
nucleon mass at $ = 0.) 

In all cases shown here, the nucleon pole exists only up to some value 
of the scalar potential. Fig.l shows that, in comparison to the exact Fad- 
deev calculation, the result of the static approximation decreases much too 
fast with increasing scalar potential. In particular, since the quark-diquark 
interaction in the static approximation is proportional to 1/M, it develops 
an unphysical singularity as M ^ 0. The upper dashed line in Fig.l shows 
the result if we fix the strength of the quark-diquark interaction as I/Mq. 
This agrees much better with the exact Faddeev result for small $, but for 
larger $ the attraction is underestimated. Since the main purpose of this 
paper is to present a more qualitative rather than quantitative discussion 
of the mechanism which leads to saturation in the NJL model, we will use 
a simple interpolation between these two extreme cases (1/M and I/Mq) 
for the quark-diquark interaction, which avoids the unphysical singularity as 
M — > and reproduces the exact Faddeev results for the case of the Eu- 
clidean sharp cut-off quite welQ The middle dashed line in Fig.l shows the 
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One could also introduce a momentum dependence into the quark-diquark interaction 
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result obtained by using the function — --— — with c = 700MeV to in- 

^ ^ Mo M + c 

terpolate between the two extreme cases. In the following we will frequently 
refer to this interpolation (denoted case (ii)), but the cases where the quark- 
diquark interaction is given by I/Mq (case (i)) or 1/M (case (iii)) will also 
be discussed. 

In Fig. 2 we show the resulting binding energy per nucleon as a function 
of the density for the three cases (i), (ii), (iii) discussed above. Irrespective 
of the choice for G^^ there is no stable state for cases (i) and (ii). The 
behavior of these curves is actually very similar to the one for the case of 
elementary nucleons in the sense that their curvature is negative and 
there is no saturation in the range of densities where solutions exist. For 
case (iii) some kind of saturation occurs, but it involves an unreasonably 
small nucleon mass (see Fig. 3) and large binding energy and density. This 
case resembles more the situation of the abnormal state found for elementary 
nucleons As is clear from the discussion in sect. 2. 3, a negative curvature of 
the function Mn{M), such as that shown in Fig.l for case (iii), corresponds 
to an attractive contact interaction in addition to the attractive vacuum 
contributions, and this reinforces the trend towards instability of the normal 
state and a consequent transition to the abnormal state. Moreover, this case 
shows a pathological behavior for small densities which is caused by the fact 
that with this choice of parameters (large G^) the zero density value of /o,l=o 
is positive (i.e., dominated by the u meson exchange term), see Eq. (|2.29| ). 

Figure 3 shows the effective quark mass M* (dashed lines) and nucleon 

by some kind of form factor. However, also in this case one has to fit the parameters 
of the form factor to the exact Faddeev resuh, leading to a similar situation as for our 
interpolated quark-diquark interaction strength. 
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mass (solid lines) corresponding to the three cases (i), (ii) and (iii), and 
in Fig.4 we plot the "effective potential" {S{M) — M^op) as a function of M 
for the cases (ii) and (iii) for three values of the density. For case (ii) there is 
no minimum in the effective potential for densities beyond ~ 0.22/m~^ in the 
range of M values where solutions for the nucleon mass exist, and in case (iii) 
the minimum rapidly moves towards small M with increasing density. These 
curves show the de-stabilizing effect of the attractive vacuum contributions 
( |2.51| ) to the effective a mass : For larger scalar potential (smaller M) the 



curvature of the lines in Fig. 3 decreases, the minima become flatter and 
eventually disappear. 

As we will see later, essentially the same situation is found with the 
proper-time cut-off scheme, and Appendix A shows a similar situation for 
the sharp 3-momentum cut-off scheme. For the case of elementary nucle- 
ons, it has been pointed out in several papers ^ that saturation can be 
achieved by including higher order interaction terms, which eventually make 
the coupling constants appearing in the original model density dependent, 
leading to repulsive effects. Actually, in Appendix A we show that in the 
present NJL model calculation the use of an 8-fermi interaction of the form 



Cx — Gx 



y -{ip^,T^pf]{ipYi^y (3.1) 



at the mean field level effectively leads to a renormalized vector coupling 
which increases with increasing scalar potential such that saturation can be 
achieved. Without sound guiding principles, however, the introduction of 
higher order interaction terms is somewhat unsatisfactory and will not be 
further discussed here0. 
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We have also checked that the problem of saturation cannot be solved by introducing 
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Qualitatively the same instability has already been observed in the case of 
the a or NJL model for elementary fermions Actually, we can expect 

from Fig.l that we cannot improve the situation in the present NJL model 
calculation for nucleons with internal quark structure: As we have pointed 
out in sect. 2. 3, a qualitative difference (stabilization of the system) compared 
to the case of elementary nucleons could emerge if the curve Mn{M) became 
flat for larger scalar potentials. However, Fig.l demonstrates that as long as 
there exists a quark-diquark threshold (or a quark-quark threshold at 3M if 
the diquark is unbound), there is only little room for this effect to work0. 

The presence of unphysical thresholds is a consequence of the lack of 
confinement in the NJL model. In particular, the thresholds force the nucleon 
mass to approach zero as the constituent quark mass goes to zero. In contrast 
to this situation, in models with confinement (bag models etc.) the nucleon 
mass is nonzero even (if the mass of the constituents is zero) because of the 
finite kinetic energy of the quarks inside the cavity. From this viewpoint, 
the presence of an appreciably positive scalar polarizability of the nucleon at 
larger scalar potentials can be considered as a manifestation of confinement. 

We will therefore consider a simple method to avoid unphysical thresh- 
olds, which was proposed in Ref . . This consists of introducing an infrared 
(IR) cut-off (yu) in addition to the UV cut-off (A) in the proper time regular- 
ization scheme. (Some formulae are collected in Appendix B, where it is also 

a kind of glueball field such as to make the NJL Lagrangian scale invariant , nor by 
introducing the t'Hooft 6-fermi (determinant) interaction due to a nonzero strange quark 
condensate Q. 

^•^It is possible to find cases where Mjv(M) has a positive curvature (e.g; the case of 
the 3-momentum sharp cut-off with Mq ~ 350Mey), but not large enough to lead to 
stabilization. Moreover, in these cases the binding energy of the nucleon decreases and 
its size becomes larger and larger with increasing scalar potential, which is also physically 
unacceptable. 
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shown explicitly that the thresholds of the quark-quark and quark-diquark 
loops are removed by the IR cut-off.) In the following discussion of the nu- 
merical results, we will concentrate on the comparison between the cases 
fi = 0, for which we will observe essentially the same situation as for the 
Euclidean cut-off scheme discussed above, and fi > 0. Therefore we confine 
ourselves to the case (ii) discussed above - i.e., the case where the strength 
of the quark-diquark interaction to be used in the static approximation is 
adjusted to agree quite well with the exact Faddeev result in the Euclidean 
cut-off scheme[3 Very similar results are found for case (i), but not for 
case (iii) which suffers from the unphysical singularity of the quark-diquark 
interaction as discussed above. We will fix the value of the IR cut-off as 
/i = 200MeV. (The results are very similar as long as > lOOMeV.) Once 
fi is fixed, the other parameters are determined in the same way as discussed 
above, and their values are listed in Table 1. 

Fig. 5 shows the plot for Mn{M) for the two cases /i = (lower solid 
curve) and fi = 200MeV (upper solid curve). Also shown are the diquark 
mass Md{M) and the sum Miy{M) +M for the two cases. It is interesting to 
observe that, although the IR cut-off also cancels the threshold in the quark- 
quark loop, the behavior of the diquark mass Md{M) is almost unchanged 
by the IR cut-off (see Appendix B), but the behavior of the nucleon mass 
Mn{M) is changed drastically. One could interpret this by saying that the 
confining mechanism is not at work for the quarks in the diquark but only for 

^"^Unfortunately, exact Faddeev results for the proper time regularization scheme are not 
yet available. We should also note that the IR cut-off removes the unphysical thresholds 
only in the proper time regularization scheme. Besides this, the proper time scheme has 
certain advantages, such as preserving gauge invariance, but it has the disadvantage that 
in general it leads to poles of Green functions in unphysical regions of the complex plane, 
which can cause difficulties - e.g; in the calculation of form factors. 
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the system as a whole, which is physically reasonable. The case n = 200MeV 
shows an appreciable scalar polarizability of the nucleon for larger scalar 
potentials, and accordingly the nucleon mass does not tend to zero with 
decreasing constituent quark mass. As we have discussed in detail in sect. 2. 3, 
this implies a repulsive contribution to the sigma meson mass and a reduction 
of the aNN coupling in the medium, working towards stabilization of the 
system. 

In Fig. 6 we show the binding energy per nucleon for the same two cases: 
II — (dashed line) and fi — 200MeV (solid line). The corresponding values 
of Gu, (see Table 1) have been adjusted such that the curves pass through the 
point (p, Eb/A) = (0.16/m-^ -15MeV). The result for // = (dashed line) 
is very similar to the one obtained with the Euchdean sharp cut-off (the hne 
labelled case (ii) in Fig.3) or the 3-momentum cut-off (Appendix A). The 
case fi = 200MeV, on the other hand, leads to an equation of state which 
saturates. The fact that the calculation docs not reproduce the empirical 
saturation point (the saturation density is too high) is related to the fact 
that we have only one free parameter (Go;) in the finite density calculation. 
Nevertheless, the observation that the inclusion of confinement aspects leads 
to a nucleon mass which has positive curvature (scalar polarizability) as a 
function of the scalar potential, which in turn leads to a saturating binding 
energy, indicates that the long standing problem of matter stability in chiral 
models with an effective vacuum potential of the Mexican hat shape can be 
solved by taking into account the quark structure of the nucleon. 

The effective masses of the quark and the nucleon are shown in Fig. 7 
for the same two cases (/x = 0: dashed lines, = 200MeV: solid lines). 
As we will discuss in more detail later, their behavior can be understood 
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from the relations ( p.55| ) and ( p.56| ), in terms of the density dependence of 
the quark-nucleon and nucleon-nucleon interaction. For the case /i = 0, the 
attraction due to the a meson exchange increases with increasing density 
due to the tachyon pole in the a propagator, while for the case /i = 200MeV 
it decreases due to the repulsive aNN contact term and the reduced aNN 
coupling, similarly to the case of phenomenologically successful non-chiral 
models ^ 0. Fig.8 shows the "effective potential", £{M) - MmP, for 
the two cases as functions of M for three values of the density. We see that 
for any fixed density the curvature of the effective potential decreases much 
faster with increasing scalar potential for the case /i = than for the case 
^ > 0. This demonstrates again the important stabilizing role played by the 
aNN contact term. 

Let us now analyze the difference between the cases /i = and fi > 
in more detail. In Fig. 9 we show the density dependent contributions to 
the squared a mass (Eq. (|2.52|) ) for p = 0.16/m~^ as functions of M. The 
solid lines show the "Z-graph" contributions multiplied by a factor 10, and 
the dashed lines refer to the contact terms. By comparing the two solid 
lines we clearly see the suppression of the Z-graph contribution caused by 
the reduced slope of Mn{M) (reduced coupling constant g^) for the case 
/i > 0. The largest and most important contribution, however, is due to 
the contact term in the case /i > 0, which increases with increasing scalar 
potential due to the increasing curvature of Mi^{M). The role of this term is 
further illustrated in Fig. 10, which shows the a propagator in the medium at 
zero momentum ((Jrp-) , dashed lines) and the Landau-Migdal parameter 



/i=o (Eq.( p.54D , sohd lines), for p = 0.16 fm~^ as functions of M. Since for 
the case p = the density dependent contributions to the a mass are rather 
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small, the vacuum term plays the dominant role, i.e., there is a tachyon 
pole in the a propagator, and the attraction due to a meson exchange grows 
without limits with increasing scalar potential, causing the collapse of the 
system. In the case yU > 0, however, the tachyon pole is avoided because of 
the large repulsive contribution associated with the contact term, leading to 
a smooth decrease of the a mass with increasing scalar potential. Moreover, 
due to the reduction of the aNN coupling constant, the a meson exchange 
contribution to the Landau-Migdal parameter /l=o now decreases smoothly 
with increasing scalar potential, similarly to the case of non-chiral models. 

Fig. 11 shows the same quantities at M = M*(p) as functions of the den- 
sity. For /i > the a mass in the medium is almost constant, and /£=o 
changes smoothly from attraction at small densities (dominance of the a me- 
son exchange term in Eq.( |2.54] )) to repulsion at larger densities (dominance 
of the uj meson exchange term). It is important to realize that the density 
dependence of /2=o directly related to the saturation of the binding energy 
per nucleon {Eb/A = £* / p — M^o)- Since for small densities the curvature 
of Eb/A has to be negative (see Eq. (|2.29|) ), the necessary condition for satu- 
ration is that the curvature of Eb/A turns from negative to positive at some 
density pc, which is the well-know liquid-gas phase transition. From this it is 
easy to show|^ that the necessary conditions for saturation can be expressed 
as /2^q(pc) < and (-"j^) > 0. It is clear from Fig. 11 that the case 
p = does not satisfy this condition. 

Finally, in Fig. 12 we plot the one-body cr-term, S*^, and the two-body 

^'''To see this, one expands both side of the relation p[EbIA)" — ~2{Eb/A)' + £" around 
p = Pc up to first order in (p — pc), and uses £" = n'^Vp/2pp + /£^o as well as the fact 
that at p = Pc the function {Eb /A)' has a minimum. (Here the primes indicate derivatives 
w.r.t. the density.) 



34 



cr-term, S*^^, defined in ( |2.31| ) by the solid and dashed hnes, respectively, 



and also the sum \r + ^Zmn "which appears in eqs. ( 2.42 )-( 2.44 ) by 
^Pjr 

the dot-dashed lines as functions of the density for the cases fi = and 
fi > 0. Referring to Eg. ( f2.45[ ), we understand that the decrease of S*^ with 
increasing density for the case // > is caused by the decrease of the aNN 
coupling constant, while the sharp increase for the case /i = reflects the 
increasing attraction due to the decreasing a mass (tachyon pole). The zero 
density values in both cases are S^^^.o — 2.7 x m x 0.7 ~ 32MeV", where 
the factors 2.7 and 0.7 correspond to (w)m-a/ ^ Eq. (^M) , 

respectively. The fact that this is too small compared with the experimental 
value Stttv.o = 45 ± 5MeV is partially related to the fact that in the presently 
used proper-time regularization scheme we are in the "firmly broken" regime 



3^ where < 1, and partially to the fact that we did not include the 
effect of the pion cloud The sign of the two-body a term (md/2^Q/dm) 
is determined by the sign of (5/^=0 /f^^)M=M*- -^^^ small densities, /l=o(^) 
exhibits the singularity due to the tachyon pole; i.e., it decreases rapidly 
with decreasing M. For the case /i = 0, this behavior also persists for larger 
densities (see Fig. 10), and therefore (c?/l=o/<9M)^^^^, is positive, increasing 
rapidly with increasing density. For the case /i > 0, however, the large contact 
term removes the tachionic behaviour (see Fig. 10), and {dfL=o/dM)j^^^j^^.j^, 
becomes a smooth function changing sign at p ~ 0.07/m~^. This is reflected 
by the behavior of S*^^ shown in Fig. 12. 
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4 Summary and conclusions 



In this paper wc investigated whether the old problem of matter stability in 
chiral field theory models, based on a linear realization of the symmetry, can 
be solved by taking into account the quark structure of the nucleon. The 
basic problem in these models is associated with the Mexican hat shape of 
the vacuum effective potential, which implies strongly attractive contribu- 
tions to the a meson mass from vacuum fluctuation effects. This decrease 
of the u meson mass with increasing scalar potential leads to an attrac- 
tive Landau-Migdal interaction between the nucleons which increases with 
increasing density. Since the necessary condition for saturation in any rel- 
ativistic mean field theory is that the attraction should decrease at high 
densities (decoupling of the a and lo mesons), it follows that stable normal 
matter cannot be described. If the nucleons are treated as elementary fields, 
it is possible to stabilize the matter by taking into account effects involving 
the polarization of the Dirac sea of the nucleon, but these give extremely large 
contributions which change the overall physical picture completely, throwing 
doubts on the reliability of this treatment. 

To incorporate the structure of the nucleon into the calculation of the 
nuclear matter equation of state, we used the Nambu-Jona-Lasinio (NJL) 
model to describe the nucleon as a quark-diquark state. We then employed 
a hybrid description of nuclear matter, in which the nucleons are moving in 
self consistent scalar and vector fields which couple to the quarks, and the 
polarization of the Dirac sea of quarks was taken into account. By consider- 
ing various expansions in powers of the density, as well as the Landau-Migdal 
interaction between the nucleons, we have demonstrated that this description 
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is very similar to that of the hnear a model for elementary nucleons. In fact, 
there is essentially only one important difference. If the scalar field acts on 
the quarks instead of directly on the nucleons, the nucleon mass need not 
be a linearly decreasing function of the scalar potential - i.e., there can be 
a curvature, which is equivalent to a "scalar polarizability" of the nucleon. 
The presence of a positive curvature of the nucleon mass as a function of 
the scalar potential is also suggested by confinement. The constituent quark 
mass decreases linearly with the scalar potential, but if the quarks are con- 
fined in the nucleon the mass of the nucleon will not decrease towards zero 
simultaneously with the quark mass. If such a scalar polarizability of the 
nucleon exists, it has two important consequences. First, there will be a re- 
duction of the (tNN coupling constant, and second a repulsive contribution 
to the cr meson mass. Both effects will increase with increasing scalar po- 
tential (density), and work against the attractive vacuum fluctuation effects 
discussed above - i.e., towards the stabilization of the system. 

One of the shortcomings of usual treatments based on the NJL model 
is, however, the absence of confinement, leading to unphysical thresholds for 
bound states like the nucleon. We have seen that, as long as these thresholds 
are present, little room is left for the effect discussed above, and one encoun- 
ters essentially the same situation as in the elementary nucleon case. The 
Landau-Migdal interaction becomes more and more attractive as the density 
increases, causing the collapse of the system. Therefore, in order that the 
stabilizing effect of the scalar polarizability can come into play, one has to 
avoid the unphysical thresholds. Several methods have been proposed ear- 
lier for this purpose, and for the qualitative purpose of this work we have 
chosen the most simple one, which consists in introducing an infrared cut-off 



37 



in addition to the ultraviolet cut-off, in the framework of the proper-time 
regularization scheme. We have shown that for physically reasonable values 
of this infrared cut-off the nucleon mass as a function of the scalar poten- 
tial does indeed have an appreciable curvature, in particular for large scalar 
potentials. We have demonstrated that this effect can lead to a saturating 
nuclear matter equation of state. That is, we have confirmed that effects 
based on the quark structure of the nucleon work towards a solution of the 
long standing problem of matter stability in chiral models based on the linear 
realization of the symmetry. 

Before drawing firm conclusions, however, it should be noted that for 
technical reasons our treatment is still afflicted with an ambiguity. That is, 
in order to simplify the Faddeev treatment of the quark-diquark state in the 
scalar potential, we used the static approximation to the Faddeev equation, 
where the momentum dependence of the quark exchange kernel is neglected. 
While this approximation works reasonably well in free space, in the medium 
it becomes worse due to the reduced constituent quark mass. Therefore we 
used an interpolated strength for the static quark exchange kernel, adjusted 
so as to reproduce the exact Faddeev result, which at present is available 
only for the Euchdean sharp cut-off scheme. We then used this interpolated 
form of the quark exchange kernel for other cut-off schemes as well. To 
resolve this ambiguity, one should either use the exact Faddeev equation 
or introduce an approximation which is valid for small constituent quark 
masses. Nevertheless, the results of this paper strongly suggest that in a more 
complete treatment one should be able to describe stable nuclear matter in 
a relativistic mean field treatment of chiral theories by letting the mean field 
couple to the quarks, provided that a kind of confinement mechanism ensures 
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that the nucleon mass does not tend to zero with increasing scalar potentiaL 
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Appendices 



A Results for the 3-momentum sharp cut-off 
scheme and stabihzation due to higher or- 
der interaction terms 

In this appendix we show the results obtained in the 3-momentum sharp 
cut-off scheme, including the stabilizing effect due to the 8-fermi interaction 
term ( |3.1| ). (The stabilization due to the 8-fermi interaction is, however, 
independent of the regularization scheme.) 

Defining for the purpose of this appendix the field o = {p\i\)'}\)\p) in addi- 
tion to the vector field lo^ = {p\ilj'-f^tp\p) , the inclusion of the 8-fermi inter- 



action (|3.1|) on the mean field level modifies the Lagrangian (|2.2|) to 



C = ^{t^-M- 2G^7^u;^) ^ - G^a^ + G^Uf^u^ - SGga'cu^cu^ + Ci, (A.l) 

where M = m — 2(5^0", G^^ = G^ — G^a"^ and = + G^uj^uj^^. Since the 
nucleon energy has the form ( p.l4| ) with G^i — > G^)-, the expression for the 



energy density in the hybrid model becomes (cf. Eq. (p.9|) ) 

£ = £vq{M) + G^a^ - G^uo^uj^" + ?>G^a^uo^uo^ + QG^uo^p + £f. (A.2) 

where £vq{M) is the quark loop term (the first term in Eq.( p.lO|) ), and £f is 
given by ( p.l9| ), where now kj^f = k — GGt^oj. To calculate the variation w.r.t. 
a and uj^, one has to take into account that now M, which appears in the 
quark loop as well as the nucleon Fermi motion term due to = Mn{M), 
as well as fcjv depend on both the scalar and vector fields. We obtain 
d£ d£ 

-— = -AG^uj>^a—--2G^uj^^ + QG,a''u>^ + QG^fs = ^ (A.3) 
ooj^ oM 
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d£ ~ dS 

= 2G^-^-2G^a-6G8aLu^cuf' + 12GsaLu^fB = 0, (A.4) 



where the baryon current is as defined below Eq. (|2.14|) , and the derivative 



w.r.t. M acts on the quark loop and nucleon Fermi motion terms {Svq and 
Sp). One can now show that u'^ = is still a solution. If we directly insert 
this relation into ( |A.3| ) and ( [A.4|) , both equations give the same requirement, 



a = dS/dM. Using uj'^ = to eliminate the vector field, and setting 
= for nuclear matter at rest, we can write the expression for the energy 
density in the same form as Eq.(2.17) for j'^ = 0, where in the denomi- 



,2^ 



nator of the second term in ( |2.1(J| ) should be replaced hj G-,, = G^ {1 + Pp 

with (3 = 9Gs/Gt^. To discuss the stabilizing effect of Gs, it is more instruc- 

(M — m)^ / 1 1 \ 

tive, however, to combine the shift -= with the the term 

4 \g^ gJ 

9G^p^ of ( p.l9| ). This leads to the final result that the energy density is given 



by the same expression as before, but with £^ replaced by 

4 = 9G.p^ (l - ^^^^^^) , {P = 9Gs/G.). (A.5) 

(Mq — m)^ Gs 

If we impose the condition 1 — — — > 0, then the term in the 

4G^ Ga, 



brackets of ( |A.5| ) is always positive for M < Mq. (This condition corresponds 



to G^ > 0.) The important point is that the repulsion ( |A.5| ) increases with 
decreasing M, or in other words, with increasing scalar potential (density). 
If Gg is chosen large enough, this term will ensure that the curves in Fig. 4 
will go up on the low M side and the minimum in the effective potential will 
not disappear. The calculation shows that this is the case for (3 > Ifm^. 

In Fig. 13 we show the nucleon and diquark masses as functions of the 
scalar potential for the 3-dimensional sharp cut-off scheme, using case (ii) 
for the interpolated quark-diquark interaction strength as discussed in sect. 
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3. (The parameters are given in Table 1.) The behavior of the sohd hne 
in Fig. 13 is very similar to the one labeled by (ii) in Fig.l and therefore we 
cannot obtain a stable nuclear matter ground state for Gs = 0, as shown by 
the dashed line in Fig.l4. The solid line shows the result for (3 = 1.5/m^. 
We wee that the binding energy indeed saturates, but at a too low density. 
(It turns out that this cannot be improved by choosing a different value for 
p.) In Fig. 15 we show the effective quark and nucleon masses for the cases 
/3 = and f3 = 1.5/m^. The stabilizing effect of the 8-fermi interaction term 
is clearly seen. 

B Some explicit formulae for the proper time 
regularization scheme 

To apply the proper time regularization scheme to a momentum loop integral 
one first follows the usual procedure of covariant loop integration, i.e., 
one combines the denominators by introducing Feynman parameters, intro- 
duces a shift of the loop momentum and performs the Wick rotation. In the 
resulting expressions one then makes the replacements 

InA -> - / — e 

J1/A2 r 

drr^-^e""^ {n > 1), 



A" {n - 1)! A/A2 

where A depends on the momenta and Feynman parameters, and A and n 
denote the UV and IR cut-off, respectively. Applying this procedure to the 
quark loop term in the energy density (first term in (|2.10|) ), the scalar qq 



bubble graph (|2.5|) and the quark-diquark bubble graph (|2.8|) , we obtain the 
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expressions 



£vq = ^2 {CsiM') - C,{M',)) (B.l) 
= -^(2C2{M') + q'Ci{M')-q' f'dxx{l-2x) 

iOTT^ V JO 

^ MMKWi ) 



X 



[mI - - g2 (1 - 2x) 



2 



+ Mj Ci(M|,) + J^dxx (^^^+M 

^A{M^,Ml)/A^ _ g-A(Af2,M|,)//i2- 



Here we defined Cnirn^) = ^e^™"^ and the and x-dependent func- 
tions A{m1,m2) = mf + (mg — mf)a; — q'^x{l — x)). The diquark t-matrix 
Ts has been approximated by the constant + pole terms, {Ts{k) ~ AiGs + 
Qs/ik"^ — M|,)), as has been explained in sect. 3. 

It is clear from the above expressions that the imaginary parts of Ug and 
Hat are canceled if /i^ > 0. It is also easy to verify that Ils{q^) < is a 
monotonically decreasing (increasing) function of (M^). Therefore, the 
condition that in (p.4| ) there is a diquark pole (M^, > 0), even for M = 0, 
is n,(g2 = 0, M2 = 0) > -1/(2G,), which in turn requires - 7r^/{3Gs) < 
^2 ^ j {^Q^y (The second inequality must be satisfied in order to have 

a non-trivial solution of the gap equation at zero density.) The parameters 
used in sect. 3 do not satisfy the first inequality, which means that the aspect 
of "confinement" (i.e., to have a finite mass of the composite system even if 
the mass of the constituents tends to zero) is not present for the diquark. It 
is, however, present for the nucleon as a whole, as is indicated by the results 
shown in Fig. 5. 
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Eucl. 


3-mom. 


pi 

= 


oper time 
^ = 200 MeV 


m[MeV] 
A [MeV] 
[GeV-2] 

Ts — Gs/G-ir 


8.99 
739.0 
10.42 
0.65 


5.96 
592.7 
6.92 
0.73 


17.08 
636.7 
19.76 
0.51 


16.93 
638.5 
19.60 
0.49 



Table 1: Parameters used for the Euclidean sharp cut-off ("Eucl.") the 3- 
momentum sharp cut-off ("3-mom.") and the proper time regularization 
schemes. For the proper time regularization scheme, we give two parameter 
sets with different IR cut offs //. 
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Figure Captions 



1. The solid and dashed hnes show the nucleoli mass Mj^ as functions 
of the 'scalar potential' $ = Mq — M in the Euclidean sharp cut-off 
scheme. The solid line is the exact Faddeev result (r^ = 0.655), and 
the dashed lines refer to the static approximation (r^ = 0.646, other 
parameter given in Table 1) for the three cases of the quark-diquark 
interaction strength discussed in the main text. Upper dashed line: 
case (i), middle dashed line: case (ii), lower dashed line: case (iii). The 
dot-dashed line shows the diquark mass Mo for the case = 0.646, 
and the double-dot dashed line shows the quark-diquark threshold for 
this case. The dotted hne shows the quark mass M — Mq — $ for 
orientation. 

2. The binding energy per nuclcon as function of the density in the Eu- 
clidean sharp cut-off scheme for the three cases of the quark-diquark 
interaction strength discussed in the main text. For each case, has 
been adjusted such that the curve passes through the point (p, Eb /A) — 
(0.16/m-^-15Mey), leading to 9G^/G^^ 5.68, 6.91 and 11.97 for 
cases (i), (ii) and (iii), respectively. 

3. The effective quark mass (dashed hnes) and nucleon mass (solid hnes) 
as functions of the density in the Euclidean sharp cut-off scheme for 
the three cases of the quark-diquark interaction strength discussed in 
the main text. Upper lines: case (i), middle lines: case (ii), lower lines: 
case (iii). The values of G^^ used in these three cases are given in the 
caption to Fig. 2. 
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4. The 'effective potential' S—M^oP as function of M for four values of the 
density in the Euclidean sharp cut-off scheme for the cases (ii) and (iii) 
of the quark-diquark interaction strength discussed in the main text. 
The values of used in these two cases are given in the capture to 
Fig. 2. The values of the density are: p — (dotted hne), p — 0.08/m~^ 
(dashed lines), p = 0.16/m~^ (sohd hues), p = 0.24/m~^ (double-dot 
dashed lines). 

5. The solid hues show the nucleon mass Mn as function of the 'scalar 
potential' (f) = Mq — M in the proper-time regularization scheme, re- 
ferring to case (ii) for the quark-diquark interaction strength. Lower 
solid line: case (a) [p = 0, = 0.492), upper solid line: case (b) 
{p = 200Me\^, Vs — 0.508). Also shown are the diquark mass Md and 
the sum Md + M for these two cases. The dotted line shows the quark 
mass M — Mq — $ for orientation. 

6. Binding energy per nucleon as function of the density in the proper- 
time regularization scheme, referring to case (ii) for the quark-diquark 
interaction strength. Shown are the results for the following two cases: 
(a) = (r, = 0.492, 9G^/G^^ 3.971, dashed line), (b) p = 200Mey 
(r^ = 0.508, 9G^/G^^ 3.323, solid line). 

7. Effective quark and nucleon masses as functions of the density for the 
cases (a) (sohd lines) and (b) (dashed lines) described in the caption 
to Fig.6. 

8. The 'effective potential' £ — M^oP as function of M for four values of 
the density in the proper-time regularization scheme, referring to the 
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cases (a) and (b) described in the caption to Fig. 6. The values of the 
density are as in Fig. 4: p — (dotted line), p — 0.08/m~^ (dashed 
hnes), p — (solid lines), p — 0.24: fm"^ (dash-double dotted 

hues). (The curve for the case p = is practically the same for the 
cases (a) and (b), and the dotted hue here refers to case (b).) 

9. The solid hnes show the "Z-graph" contributions to multiphed by 
a factor 10, and the dashed lines show the contributions of the contact 
term to as functions of M for the density p — 0.16/m~^. The cases 
(a) and (b) refer to the two cases described in the caption to Fig. 6. 

10. The dashed hnes show and the sohd lines show /i=o as func- 

tions of M for the density p — 0.16fm~^. The cases (a) and (b) refer 
to the two cases described in the caption to Fig. 6. 



functions of the density. The cases (a) and (b) refer to the two cases 
described in the caption to Fig.6. 

12. The sohd lines show S*jv, the dashed lines show S*jvAr) the dash- 
dotted hnes show S*jy^-|-^E*jv as functions of the density. The cases 
(a) and (b) refer to the two cases described in the caption to Fig.6. 

13. The solid hne shows the nucleon mass Mjv as function of the 'scalar 
potential' $ = Mq — M in the 3-momentum sharp cut-off scheme, 
employing the static approximation to the quark exchange kernel and 
the strength of the quark-diquark interaction as in 'case (ii)' described 
in sect. 3. The dashed line shows the diquark mass Mp, and the dot- 




M=M* 



, and the sohd lines show /£^o 
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dashed line shows the quark-diquark threshold. The dotted line shows 
the quark mass M — Mq — <l> for orientation. 

14. The binding energy per nucleon as function of the density in the 3- 
momentum sharp cut-off scheme, using the nucleon mass as shown by 
the solid line in Fig. 13. The dashed line shows the result for Gs — 
0, and the sohd hne refers to the case QCg/G^ = 1.5/m^. For each 
case, has been adjusted such that the curves passes through the 
point {p,Eb/A) = (0.16/m-^-15Mey), leading to 9G^/G^^ 9.84 
and 22.01 for the cases Gs — and dGs/G^^ = 1.5/m^, respectively. 

15. The effective quark mass (dashed hnes) and nucleon mass (solid hnes) 
as functions of the density in the 3-momentum sharp cut-off scheme for 
the cases Gs — and dGg/G^^ — 1.5/m^. 
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